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8  Abstract 

9  Strong  surface  waves  and  currents  generated  by  major  hurricanes  can  produce 

10  extreme  forces  at  the  seabed  that  scour  the  seafloor  and  cause  massive  underwater 

1 1  mudslides.  Our  understanding  of  these  forces  is  poor  due  to  lack  of  concurrent 

12  measurements  of  waves  and  currents  under  these  storms.  Using  unique  observations 

13  collected  during  the  passage  of  a  category-4  hurricane,  Ivan,  bottom  stress  due  to  currents 

14  and  waves  over  the  outer  continental  shelf  in  the  Gulf  of  Mexico  was  examined.  During 

15  the  passage  of  Ivan,  the  bottom  stress  was  highly  correlated  with  the  wind  with  a 

16  maximum  of  about  40%  of  the  wind  stress.  The  bottom  stress  was  dominated  by  the 

17  wave-induced  stresses,  and  exceeded  critical  levels  at  depths  as  large  as  90  m. 

18  Surprisingly,  the  bottom  damaging  stress  persisted  after  the  passage  of  Ivan  for  about  a 

19  week,  and  was  modulated  by  near- inertial  waves. 
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20  1.  Introduction 

21  Hurricanes  can  produce  extreme  forces  at  the  ocean  bottom,  even  on  the  outer 

22  continental  shelf.  Bottom  stresses  resulting  from  near  bottom  flows  are  less  noticeable 

23  than  surface  winds  and  waves.  These  episodic  wind  events  modify  and  control  the  near- 

24  bottom  environment  through  resuspension  and  transport  of  sediments,  and  redistribution 

25  of  organisms  and  chemicals.  In  the  presence  of  strong  surface  waves,  the  combined 

26  current-wave  stress  is  considerably  larger  than  the  bottom  stress  associated  only  with 

27  mean  currents  [e.g.,  Grant  et  al.,  1984;  Madsen  et  al.,  1993;  Cacchione  and  Drake, 

28  1982].  Stresses  generated  by  the  surface  waves  over  the  sea  floor  of  the  continental  shelf 

29  have  been  often  underestimated  since  wave  heights  under  hurricanes  were  believed  to 

30  exceed  20  m  or  so  only  in  the  100-year  storm  events.  Recently,  an  extreme  wave  with  a 

3 1  crest-to-trough  height  of  28  m  was  measured  under  Hurricane  Ivan  and  was  not 

32  considered  a  rogue  wave  height  but  rather  a  common  wave  height  that  can  occur  under  a 

33  major  hurricane  [  Wang  et  al.,  2005].  Significant  wave  heights  likely  surpassed  21  m  and 

34  maximum  individual  wave  heights  may  have  exceeded  40  m  near  the  eyewall  [Wang  et 

35  al,  2005]. 

36  The  Gulf  of  Mexico  (GOM)  region  provides  nearly  30%  of  the  United  States  oil 

37  supply  and  20%  of  its  natural  gas.  Hurricanes  are  major  threats  to  the  integrity  of  offshore 

38  operations  over  the  GOM  outer  continental  shelf.  Significant  damage  can  occur  to 

39  underwater  pipelines  and  to  other  underwater  infrastructures  such  as  oil  and  gas  platforms 

40  [Cruz  and  Krausmann,  2008],  There  are  reportedly  at  least  50,000  km  of  pipeline  on  the 

41  seafloor  of  the  GOM  [MMS,  2006].  Damage  to  pipelines,  which  often  is  difficult  to 

42  detect  unless  the  damage  is  catastrophic,  can  be  more  costly  to  repair  than  damage  to  the 
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43  superstructures  on  platforms.  Major  oil  leaks  from  damaged  pipelines  could  have 

44  irreversible  impacts  to  the  ocean  environment.  Improved  understanding  and  accurate 

45  prediction  of  hurricane-induced  bottom  stresses  on  the  seafloor  area  along  hurricane 

46  paths  can  enable  better  engineering  designs  to  reduce  pipeline  failures.  However, 

47  understanding  of  hurricane-induced  extreme  bottom  stresses  is  hampered  by  the  lack  of 

48  direct  measurements  of  near-bottom  flow  generated  by  winds  and  waves  under  intense 

49  storms.  Deployments  of  large  numbers  of  instruments  along  hurricane  tracks  are  not 

50  practically  feasible.  There  are  only  a  few  reported  wave  and  current  measurements 

5 1  directly  under  the  paths  of  historical  hurricanes.  For  example,  elevated  bottom  stress  and 

52  sediment  resuspension  over  the  Mid-Atlantic-Bight  shelf  off  the  east  coast  of  US  were 

53  found  during  the  passage  of  Hurricane  Edouard  in  1996  [Dickey  et  al.,  1998;  Chang  et 

54  al.,  2001].  The  eye  of  the  hurricane  passed  within  110  km  of  the  mooring  site  when  the 

55  maximum  wind  speed  and  wave  height  were  about  20  m  s'1  and  9  m,  respectively,  and  the 

56  resuspension  of  sediments  was  up  to  30  m  above  the  seabed.  The  maximum  bottom  stress 

57  based  on  a  current- wave  interaction  model  [Chris toffersen  and  Jonsson,  1985]  was  about 

58  0.35  N  nT2. 

59  Measurements  of  bottom  pressures  and  full  water-column  current  profiles  were 

60  fortuitously  made  under  the  eye  of  Hurricane  Ivan  in  2004  by  the  Naval  Research 

61  Laboratory  (NRL)  as  part  of  the  Slope  to  Shelf  Energetics  and  Exchange  Dynamics 

62  (SEED)  project  [Teague  et  al.,  2007].  Six  current-profiler  moorings  along  with  wave-tide 

63  gauges  were  deployed  on  the  continental  shelf  at  water  depths  ranging  between  60  and  90 

64  m  near  29.4°N,  88°W  (Figure  1).  Ivan  passed  over  the  SEED  mooring  array  on 

65  September  16  around  0000  UTC  as  a  category-4  storm  before  making  landfall  near  Gulf 
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66  Shores,  Alabama  [Teague  et  al.,  2007;  Powell  et  al.,  1998].  The  measurements  indicated 

67  that  significant  wave  heights  exceeded  20  m;  near-bottom  wave-induced  oscillatory 

68  currents  were  over  2  m  s'1  [ Wang  et  al.,  2005];  and  bottom  scours  exceeded  0.3  m  at  the 

69  60  m  isobath  [Teague  et  al.,  2007].  Extreme  waves,  currents,  and  scours  as  observed 

70  during  Ivan  are  likely  produced  by  other  hurricanes. 

7 1  The  main  objective  of  this  paper  is  to  quantify  bottom  stress  over  the  continental 

72  shelf  during  the  passage  of  Ivan.  The  seabed  frictional  coefficients  and  bottom  stresses 

73  were  evaluated  from  the  observed  near-bottom  currents  and  wave-orbital  velocities  in 

74  combination  with  the  wave-current  boundary  layer  model  [e.g.,  Christoffersen  and 

75  Jonsson,  1985],  This  study  focuses  on  the  growing  and  relaxation  stages  of  the  hurricane. 

76  During  the  growing  stage  of  Ivan  (between  September  10  and  17,  2004),  the  coastal 

77  ocean  was  directly  forced  by  accelerating  winds  from  10  to  50  m  s  ,  and  during  the 

78  relaxation  stage  (between  September  18  and  25,  2004),  surface  winds  became  weak  but 

79  the  circulation  was  dominated  by  wind-generated  near-inertial  currents  [Mitchell  et  al., 

80  2005].  The  correlation  of  the  bottom  stress  with  the  surface  wind  stress  is  also 

8 1  investigated. 

82 

83  2.  Methods 

84  A  number  of  models  have  been  formulated  to  evaluate  wave-induced  bottom  stress 

85  [e.g-,  Grant  and  Madsen,  1979;  Trowbridge  and  Madsen,  1984;  Christoffersen  and 

86  Jonsson,  1985;  Glenn  and  Grant,  1987].  Some  advanced  models  [Trowbridge  and 

87  Madsen,  1984;  Glenn  and  Grant,  1987]  address  complex  seabed  sediment  conditions 

88  such  as  armoring  and  moveable  beds.  In  the  following  we  used  a  simple  fixed  flat-bed 


5 


89  model  described  in  Christoffersen  and  Jonsson  [1985]  (hereinafter  CJ85),  mainly  because 

90  the  exact  nature  of  the  seabed  sediment  conditions  were  not  known  during  the  passage  of 

91  Ivan.  CJ85  is  similar  to  the  most  commonly-used  model  [ Grant  and  Madsen,  1979],  but 

92  utilizes  an  iterative  approach  for  computing  frictional  coefficients  based  on  bottom 

93  currents,  wave  statistics,  and  sediment  grain  sizes.  CJ85  was  used  to  evaluate  bottom 

94  stress  and  sediment  resuspension  over  the  Mid-Atlantic-Bight  shelf  off  the  east  coast  of 

95  United  States  in  the  wake  of  Hurricane  Edouard  and  Hortense  in  1996  [ Dickey  et  al., 

96  1998;  Chang  et  al.,  2001].  The  selection  of  a  particular  model  would  not  make  a 

97  significant  difference  in  the  general  conclusions  presented  here.  However,  the  estimated 

98  bottom  stress  during  the  storm  should  be  treated  as  a  lower  bound,  since  the  bottom  could 

99  be  a  movable  bed,  which  tends  to  generate  a  higher-roughness  scale  than  what  was  used 

100  in  the  fixed-bed  representation  of  the  critical  bottom  stress  [e.g.,  Madsen  et  al.,  1993]. 

101  Comparison  of  several  wave-current  boundary-layer  models  are  given  in  [Souls by  et  al., 

102  1993], 

103  The  iterative  procedure  described  in  CJ85  computes  the  wave  friction  factor  ( fw  ), 

1 04  the  current  friction  factor  ( fc  ),  and  then  the  combined  current-wave  stress  (  tcw  )  on  the 

105  sea  bed,  where  rcw  =  0.5  fwp  u;Fm  ;  p  is  the  density  of  sea  water;  m  describes  the  relation 

106  between  the  current  and  the  wave,  and  is  a  function  of  fw ,  f  . ,  wave  orbital  velocity  ( uw  ), 

107  current  velocity  ( uc  ),  and  angle  between  the  current  and  the  wave;  fc  depends  on  the 

108  Nikuradse  roughness  ( kN )  [Nikuradse,  1933],  the  apparent  roughness  (kA ),  and  the  height 

109  of  the  current  boundary  layer  ( h  ).  The  bottom  stress  due  to  the  mean  current 

110  is:  rc  =  0.5  fcp  uc  .  The  combined  current-wave-dissipation  rate  ( scw  )  and  the  current- 

111  dissipation  rate  ( s c )  were  estimated  using  the  following  relationships: 
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112  scw  —  (rcw  /  pf'2  /(/<■  z)  and  sc  =  (tc  /  p)  '1  /(at  z)  ,  where  z  is  the  height  from  the  bottom, 

113  and  k  is  the  von  Karman  constant  (0.4). 

114  The  near-bottom  flow  field  was  composed  of  a  background  mesoscale  component 

115  driven  by  a  basin-scale  wind  stress  curl  in  addition  to  low-frequency  currents  ranging 

116  from  inertial  to  sub-inertial  flows  \Teague  et  al.,  2007].  The  most-significant  bottom 

1 17  currents  were  directly  wind  driven  with  superimposed  high-frequency  oscillatory  currents 

118  driven  by  surface  waves  and  swells.  Bottom  stresses  and  dissipation  rates  were  computed 

119  from  background  currents  and  wave-orbital  velocities  at  the  seafloor.  The  bottom  wave- 

120  orbital  velocities  were  estimated  from  5 12-s  burst-sampling  records  of  wave-induced 

121  dynamic  pressure  measurements  at  0.25  m  above  the  bottom  [  Wang  et  al.,  2005]  based  on 

122  linear  wave  theory.  Ocean  currents  at  0.25  m  above  the  bottom  were  approximated  by 

123  extrapolating  near-bottom  velocities  from  acoustic  Doppler  current  profilers  (ADCPs) 

124  while  assuming  the  constant  stress  layer  with  a  logarithmic  velocity  profile  (i.e.,  “law  of 

125  the  wall”)  extends  from  the  seabed  to  the  depth  of  the  first  measured  velocity  above  the 

126  ADCP.  Typically  the  thickness  of  the  wall  boundary  layer  is  about  10%  of  the  thickness 

127  of  the  bottom  mixed  layer.  During  the  passage  of  Hurricane  Ivan,  the  current  structure 

128  was  found  to  be  frictionally  driven  with  overlapping  surface  and  bottom  boundary  layers 

129  [Mitchell  et  al.,  2005]  suggesting  that  the  entire  water  column  over  the  shelf  was  either 

130  weakly  stratified  or  well  mixed.  Therefore,  a  6-  to  9-m-thick  wall  boundary  layer  is 

131  expected  near  the  vicinity  of  the  sea  floor.  The  moorings  rested  on  a  sandy  seabed  with 

132  grain  sizes  varying  from  0.06  mm  to  1  mm  with  a  median  size  of  about  0.3  mm  [ Sawyer 

133  et  al.,  2001];  therefore  kvwas  approximated  to  be  about  0.3  mm.  The  alignment  angles 

134  between  the  currents  and  wave  orbital  velocities  were  not  known,  and  therefore,  the 
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135  currents  and  waves  were  assumed  to  be  in  the  same  direction  and  the  alignment  angle 

136  was  set  to  zero  for  the  computations.  Since  wave  statistics  were  only  sampled  every  8 

137  hours,  estimates  of  tcw  and  scw  were  limited  to  3  times  per  day.  The  bottom  currents 

138  from  the  ADCP  were  sampled  at  15 -minute  intervals,  and  therefore  rc  and  sc  were 

139  computed  at  a  higher  sampling  rate  by  interpolating  8-h  estimates  of fc  into  the  ADCP 

140  sampling  rate.  By  following  Madsen  et  al.  [1993],  the  wave-orbital  velocity  was 

141  approximated  as  the  root-mean-square  amplitude  of  a  sinusoidal  wave,  where  uw  -  V2cr,  , 

142  and  <ru  is  the  standard  deviation  of  orbital-velocity  fluctuations  based  on  the  5 12-s  burst- 

143  sampling  record  of  pressure  at  0.25  m  above  the  bottom  [ Wang  et  al.,  2005].  The  surface 

144  wind  field  over  mooring  locations  was  constructed  by  combining  NDBC  buoy  winds  with 

145  Ivan  post  storm  wind  analysis  products  [  Wang  et  al.,  2005;  Powell  et  al.,  1998], 

146 

147  3.  Observations 

148  A  dramatic  increase  in  bottom  stress  and  dissipation  rate  was  found  at  all  six 

149  mooring  locations  during  the  passage  of  Ivan.  At  the  60-m  isobath,  the  current- wave 

150  stress  ( zcw  )  was  enhanced  by  two  orders  of  magnitude  (Figure  2c)  and  the  current- wave 

1 5 1  dissipation  rate  ( ecw  )  was  enhanced  two  to  three  orders  of  magnitude  (Figure  2d)  as  the 

152  wind  speed  accelerated  from  10  m  s'1  to  50  m  s'1,  during  which  wave-orbital  bottom 

153  velocity  intensified  from  0(mm  s'1)  to  0(1  m  s'1)  (Figure  2a,  b).  At  the  90-m  isobath, 

154  wave-orbital  velocities  were  relatively  weak,  background  currents  were  strong  (Figure 

155  2g),  and  the  current- wave  stress  was  a  factor  of  two  smaller  than  that  at  the  60-m  isobath 

156  (Figure  2c, h).  Peak  bottom  stresses  at  moorings  Ml  to  M6  were  0.47,  0.58,  0.84,  0.62, 

157  0.48,  and  0.34  N  m'  ,  respectively.  Some  of  the  differences  in  current- wave  stresses  at 
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158  different  mooring  locations  were  related  to  the  undersampling  of  wave-orbital  velocities 

159  at  an  8-h  interval  and  depth  dependence  of  wave  statistics.  Out  of  the  six  moorings,  M3, 

160  which  was  located  to  the  right  of  Ivan’s  path  (Figure  1),  had  the  strongest  bottom  stress, 

161  largest  wave-orbital  velocity  (Figure  2b, c)  and  the  largest  surface- wave  height  of  about 

162  28  m  [  Wang  et  al.,  2005].  There  was  an  increase  in  background  stress  levels  following 

163  the  hurricane;  the  current  stress  ( rc),  averaged  over  the  mooring  array,  prior  to  Ivan  was 

164  7xl0'3  N  m'2  [lxlO'6,  4xl0'2]  while  that  stress  after  Ivan  was  30x10°  N  in2  [20xl0‘6, 

165  27x  1 O'  ],  where  minimum  and  maximum  stresses  are  given  in  the  brackets. 

166  The  ADCP  echo  intensity  (Figure  2e,j)  reflects  the  concentration  of  particles  in  the 

167  water  column  [Deities,  1999].  Therefore  the  timing  and  the  duration  of  the  resuspension 

168  of  sediments  (during  which  tcw>  the  critical  stress)  can  be  identified  from  the  echo 

169  intensity.  The  observed  critical  stress  is  consistent  with  the  spiking  of  echo  intensity.  The 

170  resuspension  occurred  as  the  wind  speed  exceeded  ~  15  ms'1  and  lasted  for  two  days 

171  over  both  60-m  and  90-m  isobaths.  After  the  passage  of  Ivan,  the  resuspension  continued 

172  at  the  60-m  isobath,  and  was  modulated  by  near-inertial  waves  (Figure  2e),  but  it  was 

173  considerably  weaker  at  the  90-m  isobath  (Figure  2j).  The  intensity  of  the  ADCP  echo  was 

174  enhanced  up  to  25  m  above  the  instrument,  implying  that  sediments  were  resuspended  to 

175  about  25  m  above  the  seabed  (Figure  2e,j). 

176  The  relationship  between  bottom  stress  (  tcw  ,  rc )  and  wind  speed  ( Um  )  during  the 

177  growing  stage  of  the  hurricane  was  studied  by  averaging  bottom  stresses  into  appropriate 

178  bins  as  a  function  of  the  wind  speed  (Figure  3).  The  impact  of  surface  waves  on  the 

179  bottom  stress  was  not  important  for  winds  less  than  8  m  s'1.  There  was  a  rapid  build  up 

180  inr^as  Um  increased  from  10  m  s'1  to  15  m  s'1.  The  growth  of  tcw  slowed  down  for  Um  > 
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181  20  ms'1,  and  rcw  was  approximately  proportional  to  U;A .  rc  was  also  proportional  to  U;0 

182  for  most  of  the  wind  record.  On  average  rc  was  a  factor  of  4  smaller  than  rcw  based  on 

183  the  wave-orbital  velocity,  uw  =  ~J2<JU  (Figure  3).  The  current- wave  stress  can  be 

184  approximated  as  tcw  «  4xlO~4[/120  for  winds  higher  than  15  m  s'1,  which  in  turn  implies 

185  that  tcw  «  Cts  ,  where  rs  is  the  surface  wind  stress,  and  the  proportionality  constant, 

186  C  =  4  x  1 0  4  /  paiCsD ;  pmr  is  the  density  of  air,  and  CSD  is  the  surface  drag  coefficient.  For 

187  CSD  =  (1.5- 2.5)  x  1(0  (13),  C  is  ~  0.15  -  0.2.  Sensitivity  of  rcw  to  the  strength  of  the  orbital 

188  velocity  was  examined  by  choosing  the  wave-orbital  velocity  as  significant  wave-orbital 

189  velocity  ( uws  ),  where,  um  =  2 <ju  ,  and  the  wave-orbital  velocity  as  maximum  wave-orbital 

190  velocity  ( uWM  ).  It  is  clear  from  Figure  3  that  tcw  is  highly  sensitive  to  the  magnitude  of 

191  the  wave-orbital  velocity  chosen.  The  estimate  of  tcw  for  the  maximum  wave-orbital 

192  velocity,  uWM  is  about  a  factor  of  2  larger  than  that  for  the  root-mean-square 

193  velocity,  V2crn  . 

194 

195  4.  Summary  and  Discussion 

196  During  the  passage  of  Ivan,  the  bottom  stress,  energy  dissipation  rate,  and 

197  resuspension  of  sediment  were  controlled  primarily  by  dissipative  processes  induced  by 

198  surface  waves,  whereas  during  the  relaxation  stage  of  Ivan,  wave-orbital  velocities 

199  became  small  and  the  bottom  stress,  dissipation  rates,  and  resuspension  of  sediments 

200  were  detennined  by  the  observed  bottom  current.  The  strongest  stresses  occurred  to  the 

201  right  of  the  stonn  path.  Bottom  stresses  exceeded  critical  levels  on  the  outer  continental 

202  shelf  at  depths  as  large  as  90  m.  Bottom  damaging  stresses  occurred  during  the  passage  of 
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203  the  storm  and  for  about  a  week  after  storm  passage.  The  bottom  stress  associated  with  the 

204  ocean  currents  peaked  at  wind  speeds  of  about  30  m  s'1  and  then  may  actually  have 

205  decreased  as  the  wind  speed  increased,  paralleling  the  surface  drag  reduction  at  high 

206  winds  [ Jarosz  et  al.,  2007;  Powell  et  al.,  2003].  The  current-wave  induced  bottom  stress 

207  continued  to  increase  as  the  wind  speed  increased,  and  was  about  1 5%-20%  of  the 

208  surface  wind  stress,  where  uw  =  -J2au .  The  maximum  stress  based  on  the  maximum  wave- 

209  orbital  velocity  was  found  to  be  as  large  as  40%  of  the  surface  wind  stress.  On  average, 

210  bottom  stresses  induced  by  ocean  bottom-current  interactions  with  ocean  surface- wave- 

211  related  currents  were  a  factor  of  4  larger  than  the  bottom  stresses  attributed  to  ocean 

212  bottom  currents  alone. 

213  The  occurrence  of  critical  bottom  stresses  on  the  outer  shelf  must  be  considered  in 

214  the  engineering  design  of  structures  on  the  bottom  and  in  detennining  where  pipe  lines 

215  should  be  buried  and  not  just  laid  on  top  of  the  ocean  floor.  Some  climate  models  have 

216  predicted  an  increase  in  the  frequency  of  intense  hurricanes  as  the  climate  warms  [Bender 

217  et  al.,  2010].  The  cumulative  effects  of  enhanced  bottom  stresses  and  associated  transport 

218  of  large  quantities  of  sediment  along  the  shelf  edge  could  be  a  trigger  mechanism  for  a 

219  slumping  or  mass-wasting  event  at  the  shelf  break. 

220 
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301  Figure  Captions 

302  Figure  1.  Path  of  Ivan  (dashed  red  line).  The  crosses  denote  the  center  of  Ivan  for  the 

303  time  marked  in  red.  Mooring  locations  at  60-m  and  90-m  isobaths  are  marked  in  solid 

304  circles  (M1-M6).  The  yellow  triangle  is  the  National  Data  Buoy  Center  (NDBC)  buoy 

305  42040.  Contours  indicate  bathymetry  in  meters.  The  translation  speed  of  the  hurricane 

306  was  about  6  m  s'1. 

307 

308  Figure  2. Observations  at  M3  (left  panels)  and  M6  (right  panels)  mooring  locations.  Time 

309  series  of  (a,  f)  10-m  wind  speed  Um  in  m  s'1  and  radial  distance  to  Ivan’s  center  (black)  in 

310  kilometers  [Powell  et  al.,  1998].  (b,  g)  Root-mean-square  estimate  of  wave-orbital  speed, 

311  uw  =  y[2<ju (red  dots),  where  cr„  is  the  standard  deviation  of  orbital-velocity  fluctuations 

3 12  based  on  the  5 12-s  burst-sampling  record  of  pressure,  and  ADCP  currents  uc  (black  line) 
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313  at  0.25  m  above  the  bottom.  Units  are  in  m  s'1,  (c,  h)  Combined  current- wave  stress  rcw 

314  (red  dots)  and  current  stress  rc  (black  line)  in  N  in'  .  (d,  i)  Estimated  current-wave 

315  dissipation  rate  scw  (red  dots)  and  current  dissipation  ec  (black  line)  at  0.25  m  above  the 

316  bottom.  Units  are  in  W  kg'1,  (e,  j)  ADCP  echo  intensity  as  function  of  height  (H)  in 

317  meters  from  the  bottom 

318 

3 1 9  Figure  3.  Mean  stress  at  0.25  m  above  the  bottom  vs  wind  speed  U m  during  the  growing 

320  stage  of  Ivan.  The  bottom  stress  based  on  mean  current  (uc)  is  the  solid  black  line  with 

321  diamonds.  The  combined  current-wave  stress  rcw  based  on  bottom  orbital-velocities: 

322  uw  =  42a  u ,  significant  wave-orbital  velocity  uws  =  2a  ,  and  maximum  orbital-velocity, 

323  uWM  are  marked  in  dots  (red),  squares  (blue),  and  open  circles  (magenta)  ,  respectively. 

324  Error  bars  are  given  by  the  thin  red  vertical  lines  and  denote  maximum  and  minimum 

325  values  of  zcw  for  uw  =  42a u  .  The  dashed  black  and  red  lines  are  not  mathematical  fits  to 

326  the  data,  but  represent  r  =  1 0-4 U2W  (black),  r  =  4  x  1 0  1 U,;  (red). 

327 

328 
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